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Abstract. We investigate pairing states realized at the (001) interface of a spin- 
triplet superconductor Sr2Ru04 on the basis of microscopic calculations. Because of 
a Rashba-type spin-orbit interaction induced at the interface, strong parity-mixing of 
Cooper pairs between a spin-singlet state and a spin-triplet state occurs in this system. 
There are also strong inter-band pair correlations between the spin-orbit split bands, in 
spite of the considerably large spin-orbit splitting. This is due to frustration between 
the spin-orbit interaction and pairing interactions. In this pairing state, time-reversal 
symmetry is restored, in contrast to the bulk Sr2Ru04 which is believed to be a chiral 
p + ip superconductor with broken time-reversal symmetry. It is demonstrated that, 
because of these features, the pairing state at the interface is a promising candidate 
for the recently proposed time-reversal invariant topological superconductor. 
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1. Introduction 

In unconventional superconductors, the BCS order parameter possesses internal degrees 
of freedom, which give rise to various rich physics, as explored for Helium 3 [I], Sr2Ru04 
[2l E], and heavy fermion superconductors [11 15]. For Sr2Ru04, a possible realization 
of a chiral p + ip pairing state is suggested by several experimental and theoretical 
studies [31 [HI [7]. In this p-wave pairing state, time- reversal symmetry is broken by 
orbital degrees of freedom. Because of this feature, the chiral p + ip superconductor 
bears some similarities to the quantum-Hall-effect (QHE) state [8J. For instance, in a 
chiral p + ip superconductor with open boundaries, a gapless edge mode propagating 
in only one direction appear at the boundary edges, which is in analogy with a chiral 
edge state in the QHE state. This similarity ultimately stems from the realization of 
a topological state in both of these quantum condensed phases. A topological state 
is a novel class of quantum ground state which is characterized not by conventional 
long-range order, but by a topologically nontrivial structure of the Hilbert space. In 
a topological state, there is a bulk excitation energy gap which ensures the stability 
of this state, and, as mentioned above, there are also gapless edge states which play 
an important role for transport phenomena. In the case with broken time-reversal 
symmetry such as the QHE state and the chiral p + ip superconductors, the topological 
structure is associated with the existence of a nonzero topological number, i.e. the 
Chern number [9]. Recently, another class of a topological state was theoretically 
proposed for band insulators [101 [HI [12] and experimentally observed [I3]. This 



topological state, which is called the Z2 topological insulator, possesses time-reversal 
symmetry, in contrast to the above-mentioned topological state without time-reversal 
symmetry, and is characterized by the existence of two counter-propagating gapless 
edge modes, which are associated with the Kramers doublet. These gapless edge modes 
give rise to the quantum spin Hall effect, which has been attracting recently much 
interest in connection with possible applications to spintronics. As there is similarity 
between the chiral p + ip superconductors and the QHE state, there is parallelism 
between the Z2 topological insulator and noncentrosymmetric p-wave superconductors 
[m [IS [161 [HI [HI [IH]- In these few years, many classes of noncentrosymmetric 
superconductors (NCSC), the crystal structures of which lack inversion symmetry, have 
been discovered [201 [211 l22l l23l [211 [25]. Some experimental and theoretical studies 
suggest that p-wave pairing states may be realized in certain systems of NCSC such 
as CePtsSi and Li2Pt3B [261 [23 [2H]- However, unfortunately, these NCSC are not 
suitable for the realization of the Z2 topological phase, because their superconducting 
gaps possess nodes [291 [21], from which gapless quasiparticles in the bulk appear, and 
destabilize the topological state. 

In this paper, we investigate a possible realization of the Z2 topological 
superconductivity at an interface of Sr2Ru04. We consider the (001) interface, at 
which the Rashba-type spin-orbit (SO) interaction breaking inversion symmetry may 
be induced [3^. We can also consider a setup in which a thin film of Sr2Ru04 is 
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fabricated on a substrate with a bias potential applied perpendicular to the (001) 
interface, which controls the strength of the Rashba SO interaction. Apart from the 
exploration of the Z2 topological superconductivity, such a system is interesting in that 
it is suitable for the systematic investigation on the effect of party-mixing of pairing 
states raised by broken inversion symmetry [311 132]. For the realization of substantial 
parity-mixing of Cooper pairs, the existence of attractive interactions in both spin- 
singlet and spin-triplet channels is crucially important. Sr2Ru04 is a good candidate 
for the realization of such a situation, because, according to microscopic analysis on the 
mechanism of superconductivity for this system, both the p-wave channel and the (i-wave 
channel enjoy substantially strong attractive interactions [331 [3l]. It is expected that the 
addition of the asymmetric SO interaction to this system raises the strong admixture 
of spin-singlet pairs and spin-triplet pairs. In general, the structure of the parity-mixed 
Cooper pairs is determined by competition between the asymmetric SO interaction 
and pairing interactions in each channel [311 [35l [361 137] . When the asymmetric SO 
interaction is dominant, the structure of the d- vector for the spin-triplet component 
is mainly constrained by the SO interaction to suppress pairings between two SO split 
bands, which are unfavorable when the SO split is much larger than the superconducting 
gap. However, in the case that the pairing interaction that is not compatible with the 
symmetry of the asymmetric SO interaction is dominant, the direction of the d-vector 
does not minimize the energy cost due to the asymmetric SO interaction, yielding inter- 
band pairings between the SO split bands. For the case of Sr2Ru04, the dominant 
pairing interaction exists in the p-wave channel, and the d-vector is perpendicular to 
the xy-^Aane because of the bulk SO interaction of the d-electron orbitals [31 [331 [381 [39]. 
In this paper, we consider the case that the asymmetric SO interaction at the interface 
is stronger than the bulk SO interaction, and examine how the chiral p + ip state in the 
bulk is affected by the asymmetric SO interaction, and what pairing symmetry is most 
stabilized at the interface. This is another purpose of the current paper. 

Our main results are as follows. As the asymmetric SO interaction becomes strong, 
the d-vector of the p-wave pairing is oriented to directions parallel to the xy-pla.ne. 
However, there exists strong frustration between the asymmetric SO interaction and 
the p-wave pairing interaction because of an additional anisotropic structure of the 
pairing interaction. Thus, the direction of the d-vector does not fully optimize the 
asymmetric SO interaction, inducing substantial amount of Cooper pairs between the 
two SO split bands, in spite of the SO splitting much larger than the superconducting 
gap. As a result, the stable pairing state possesses p + d wave symmetry, rather than 
s + p wave symmetry or d + f wave symmetry which is expected to be stabilized for 
the Rashba superconductors when the inter-band pairs are suppressed [311 [351 EH] . The 
notable feature of the p + d wave state is that the single-particle energy has a full gap, 
and there are no nodal excitations for small strength of the asymmetric SO interaction. 
Furthermore, since the d-vector is parallel to the xy-plane, time-reversal symmetry is 
restored, which makes sharp contrast with the chiral superconductivity realized in the 
bulk of Sr2Ru04. These two features are quite important for the realization of the Z2 
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topological superconductivity mentioned above. We can show that the pairing state 
realized at the interface is topologically equivalent to the combined state of a p + 
state and a p — ip state, which is time-reversal invariant, and supports the existence of 
counter-propagating gapless edge modes, which carry spin currents. 

The organization of this paper is as follows. In the sections 2 and 3, the pairing 
state at the interface of Sr2Ru04 is microscopically investigated on the basis of the 
scenario that the pairing interaction is caused by electron correlation effects. We analyse 
the structure of the parity-mixed pairing gap. In the section 4, exploiting the results 
obtained in the section 3, we present a possible scenario for the realization of the Z2 
topological superconductivity in this system. Discussion and summary are given in the 
last section. 

2. Model and formulation 

In this section, we introduce a low-energy effective model for superconductivity realized 
at an interface of Sr2Ru04, and present a theoretical framework used for the study on 
pairing states. 

2.1. Low-energy effective model 

In Sr2Ru04, there exist three quasi 2-dimensional orbitals 7iu4:dxy,xz,yz in RUO2 plane 
for electrons which play the most significant roles for low-energy properties. There 
are several theoretical proposals for the microscopic origin of pairing interaction in 
Sr2Ru04 [21 [3]. One of promising scenarios is that an effective pairing interaction in the 
p-wave channel is produced by higher order processes of electron-electron interaction 
through the Kohn-Luttinger-type mechanism [33j. In this scenario, which is first 
proposed by Nomura and Yamada, among the three bands, a,P and 7-band formed 
by the t2g orbitals, the 7-band originating from the dxy orbital is considered to be 
most important for the realization of the superconductivity. In this paper, we employ 
this scenario, because this appraoch enables us to calculate the transition temperature 
which is quantitatively in good agreement with experimental observations. Therefore, 
to discuss the appearance of superconductivity, it is sufficient to focus only on the 7- 
band. Although there exists the SO interaction between the t2g orbitals (we call this SO 
interaction the bulk SO interaction) which tends to direct the d-vector parallel to the 
z-axis, its effective energy scale for the pinning of the d-vector is small and negligible for 
a discussion on the transition temperature Tc. However, for electrons near the surface or 
an interface parallel to RUO2 plane, there exists another kind of spin-orbit interaction 
called the asymmetric SO interaction which breaks both reflection symmetry in the 
momentum space k ^ —k and spin rotation symmetry. This SO interaction originates 
from the spin-flip hopping processes between the Ru t2g orbitals, of which the wave 
functions are modulated by a potential gradient in the vicinity of an interface. For the 
superconductivity, it tends to make the direction of the d-vector perpendicular to the 
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2;-axis as shown in sectioE|31 In the following analysis, we assume that the asymmetric 
SO interaction is sufficiently stronger than the bulk SO interaction, and neglect effects 
of the bulk SO interaction. 

We, then, simply describe the electrons of the 7-band near an interface by the single 
band Hubbard model, 

H = ^cl[ek + aCo{k) ■ <j]ck + UY^ n^ni^, (1) 

k i 

where Ck = {ck^yCkiY is the annihilation operator and riia = c[„Cia. The asymmetric 
type SO interaction induced near the (OOl)-interface is incorporated into the first term 
of the Hamiltonian. The strength of the asymmetric SO interaction is denoted by a. 
We assume that the Rashba form of the asymmetric SO interaction [30j. For Sr2Ru04, 
the dispersion relation Sk and the Rashba type SO interaction are approximated by 

£fc = — 2ti(cOS fca; + cos /Cy) + 4it2 cos fca; cos /Cy — /i, (2) 

/^o(fc) = (sin fcj,, — sin A;^, 0). (3) 

The parameters are fixed as (^1,^2) = (1-0, —0.375) taking ti as the energy unit, and the 
filling is n = 1.32. The hopping integral and the electron density are chosen so that the 
Fermi surface of our model is consistent with the experiments. In the real system, the 
form of Colk) may be more complicated. However, as will be shown in the section 3, 
this simplified model captures important physics raised by broken inversion symmetry. 



2.2. Perturbation theory for the Kohn- hutting er mechanism of superconductivity 

There are some theoretical proposals for the mechanism of p-wave superconductivity 
realized in Sr2Ru04. One promising scenario is that the pairing interaction in this 
system is caused by the Kohn-Luttinger mechanism; higher order interaction processes 
due to the Coulomb interactions U give rise to effective pairing interactions in interaction 
channels with nonzero angular momentum [40]. Actually, Nomura and Yamada 
demonstrated that interaction processes up to the third order in U yield a strong pairing 
interaction in the p-wave channel for the microscopic model of Sr2Ru04. This scenario 
successfully explains the origin of the p-wave superconductivity realized in this system. 
We, here, apply this perturbation theory for the pairing interaction to the model ([T]). 
For this purpose, we introduce noninteracting Green's function, 

G%{k) = E l%{k)Gl{k), (4) 

T = ±l 

t^,{k)=\{l + TCo{k)-CT)^^, (5) 

Gl{k) =-^— , (6) 

ILOn — ^kr 



where Skr = ek + Ta\Co{k)\, Co{k) = Coik)/\Co{k)\ and \Coik)\ = \/Ei=i[^oiikW- 
is the fermionic Matsubara frequency. As seen in the form of G^^, the Fermi surface 
splits into two bands whose dispersions are e+ and £_ with the splitting a/vp {vp is 
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the averaged Fermi velocity of the two bands). Note that Cq = at the van Hove points 
(0, ±7r), (±7r, 0) and the Fermi surface is changed httle around them since the Rashba 
SO interaction is small there. 

The effective pairing interaction is expanded up to the third order with respect to 
U. It is expressed as, 

VaiaQasaA^, k') — 2 [Kr^CTf(T30-4 + K^CT2(T30-4 

VZ^{k, k') + UWk + k') + U%x\k - k')f + {x\k + k'))\ 

V^ZAK k') = 2U'Rej:Gl{-k + q)[x\q)Gl{-k' + q) - <P\q)Gl{k' + q)], 

VZtik, k') = - U\x\k - k') - x%k + k% 

V^ZAk, k') = 2U'ReY.GlAk + + n<l)\[GL{k' + q) - Gl^{-k' + g)], 

Vcraacr{k., k ) = V^CT(Ta(^) k), 



others = 0, 



where 



x°(?)= -^EG°.(g + ^)G'LW, (7) 

k 

-^i:GL{'l-k)Gl{k), (8) 

k 

and k = {ujn, k) , and T and are, respectively, temperature and the number of Ru sites. 
X° and 0° do not depend on spins because G'^^{k) = G^^{k) is satisfied. We, here, neglect 
many terms in V which arise from non-zero off-diagonal elements of Green's function 
Gcra{k), because G^^{k) is smaller than G^^r by the factor oia/ep where Sp is the Fermi 
energy. Besides, the terms in V with spin-fiip processes represent the perturbative effects 
of the Rashba SO interaction and, in itself, do not have crucial importance as long as 
a <^ Ep. Each ^^1^2^30-4 consists of the RPA-like terms V^^^ and the vertex-correction 
terms \/^'^''. The former is included within random phase approximation (RPA), and 
the latter is not. For spin-singlet pairing, the RPA-likc terms give dominant attractive 
interaction and play significant roles for its stability, while the vertex-correction terms 
do for triplet pairing. 

The transition temperatures for the superconductivity are calculated by solving the 
linearized Eliashberg equation 

XA^^ik) = - ^ E Vaa'pfS'iK k')Gl^{k')G%^,{-k')A,,,{k') (9) 

k 

where A is the anomalous self-energy and A is the eigenvalue. We identify the 
temperature for which A(T) = 1 as the transition temperature. The normal self-energy 
is neglected because it is not important in the present study. In this equation, spin- 
fiip processes are included only in the factor G{k)G{—k) and are not in V within our 
approximation. This is because the factor G{k)G{—k) behaves like a window function 
which allows electrons only near the Fermi surface to participate in the superconductivity 
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and therefore has non-perturbative effects of the Rashba SO interaction, while the spin- 
flip scattering processes in V are perturbative. We note that some of the elements 
of fc) are strongly anisotropic in the fc-space, which restrict the possible 

symmetries of the gap functions. Furthermore, these elements with /3 7^ or 7 7^ 7' 
give rise to parity-mixing between spin-singlet and spin-triplet pairs, which is one of the 
most remarkable features of NCSC 

The anomalous self-energy is generally written as, 

^{k) = [DQ{k) + D{k) ■ (T]ia2, (10) 

where DQ{k) and D{k) are the singlet and triplet parts, respectively, and the k- 
dependence of D^{iLJn,k) represents the symmetry of the superconductivity. The 
structure of the d-vector D is determined by the two factors. One is the pairing 
interaction V and, in the case of Rashba superconductors, there exists the other factor 
Co- The microscopic origins of these two factors are generally different, and the d- 
vector which V favors does not necessarily coincide with the one which Cq favors. For 
sufficiently large a with <^ a ^ ep, the most stable direction of the cZ-vector 
is D II jCq because, if this condition is satisfied, A(A;) in the matrix form can be 
diagonalized with respect to the r = ± bands. Conversely, when D is not parallel to Cq, 
inter-band pairs between the SO split bands are induced, which generally lead to pair- 
breaking effects. In contrast, for enough small a, the structure of D is determined by 
the pairing interaction. Generally, these two factors which can determine the structure 
of the d-vector compete with each other. In the next section, we discuss the effects of 
the Rashba SO interaction on the transition temperature and the structure of D^. 

To solve the Eliashberg equation numerically, we divide the Brillouin zone into 
64 X 64 meshes and take 1024 Matsubara frequencies. 

3. Results for the pairing state 

In this section, we present the results for stable pairing states and their transition 
temperatures calculated by using the formulation given in the section 2. 

3.1. Structure of pairing interaction 

We first show the momentum profile of figure [1] for a = and a = 0.1 at 

T = 0.01, which plays an important role for the pairing mechanism. The difference in 
for a = and a = 0.1 are so small that the fc-dependence of V is also changed little 
by the Rashba SO interaction. Indeed, we have confirmed that the fc-dependence of V 
is almost unchanged at least up to a ~ 0.1. This means that, in view of the pairing 
interaction, the most stable symmetry of the superconductivity for a = is stable also 
for a 0. As shown in figure [21 however, the amplitude of the pairing interaction 
for the triplet part Vt = ^[V^^ + V^^^„] is decreased for large a, because the main 
part of Vt is V^*^^ which is sensitive to the electronic structure. On the other hand, 
the interaction for the singlet part Vs = [V^^^ + V^^^^] is not so changed, since Vg is 
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Figure 2. Plot of the pairing interaction Vt(fc, k') on the k — {k^, ky) plane with fixed 
fc' at T = 0.0f,J7 = 5.5, for a = (left panel) and a = 0.1 (right panel). Here the 
Matsubara frequencies which appear in Vt is fixed as Vn — v'^ = inT . fc' is fixed as 
illustrated in the figure. 

mainly determined by V^^^^ which is directly related to the a-insensitive function ■ 
Therefore, it is expected that the triplet superconductivity would be suppressed while 
the singlet superconductivity unaffected through the change in V by the Rashba SO 
interaction. However, as mentioned before, the Rahsba SO interaction has the other 
important effects on which are non-perturbative in the sense that G{k)G{—k) strongly 
restricts the possible symmetries of the gap functions, and gives rise to parity-mixing of 
Cooper pairs. 

3.2. Pairing state and transition temperature in the case without parity-mixing 

As mentioned above, there are two important effects of the Rashba SO interaction on 
pairing states: one is to constrain the direction of the d-vector of spin-triplet pairings, 
and the other is parity- mixing. We, first, examine the former effect, neglecting the effect 
of parity- mixing for a while. Actually, the parity-mixing is not negligible in the present 
study, and will be discussed in the next subsection. 

Neglecting the terms which mix the singlet and the triplet gap functions in the 
Eliashberg equation ([9]), we calculated the transition temperatures for the spin- 
singlet channels and for the spin-triplet channels separately. We also computed the 
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Figure 3. a versus Tcs and Tct at U ~ 5.5. The gap functions roughly expressed 
as Dq ~ (cos fcj, — cosfcy) for T^s, D ^ cosk^sinkyX + cos ky sink^y) for T^t 
(i), D ^ [—{acoskx — bcosky) sinkyX + {bcoskx ~ a cos ky) sin k^y] for Tct ^^^d 
D ~ (cos ky sm kx + i cos k^ sin ky)z for Tct (hi). 



fc- dependence of the gap functions self-consistently from ([9]) . Figure [3] shows the a- 
dependence of Tc for the singlet (Tcs) and the triplet (Tct) superconductivity at U = 5.5. 
The solid line with closed squares is Tcs and the gap function for the spin-singlet pairing 
is roughly given by that with dx2_y2 symmetry, Dq ~ [coskx — cos ky). For the singlet 
pairing, this is the only one stable gap function, as in the case without the Rashba SO 
interaction [33]. The other lines in figure [3] are for the spin-triplet states, and calculated 
with the assumption that the d-vector belongs to (i) the Ai representation of the point 
group 04^,, (ii) Bi and (iii) E, respectively. For these representations, the d- vector is 
roughly of the form of 

(i) D^^ ~ (— cos kx sin kyX + cos ky sin kxi/), 

(ii) D^^ ~ [—(a cos k^ — b cos ky) sin kyX + {b cos k^ — a cos ky) sin k^y] with a > b, 

(iii) ~ (cos ky sin kx + i cos kx sin ky)z. 

All of these gap functions are p-wave gap functions. Note that they are all different 
from the form {coskx — cos ky)CQ for which the gap function A can be diagonalized 
with respect to the SO split bands and no inter-band pairing is realized [311 [351 1361137] . 
Thus, in the spin-triplet pairing states obtained in this calculation, there are always 
inter-band Cooper pairs. This is due to incompatibility between the symmetry of the 
pairing interactions and the symmetry of the Rashba SO interaction, as mentioned 
before. It should be also notified that the triplet states with D^^ and D^^ are time- 
reversal invariant, while the triplet state with is not, but a chiral p + ip state. 
This is easily seen as follows. For the state with D^^ or D^^ , under time-reversal 
operation, the gap function for ff pairs A^^{k) = —Di{k) + iD2{k) is transformed as 
-Di{-k) -iT>2(-fc) = A||(fc), and also, A^(fc) A||(fc). Thus, the D^^ state and 
the D^' state are time-reversal invariant. 

As can be seen in figure El Tcs is not so strongly affected by the Rashba SO 
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interaction. In contrast, is rapidly decreased, as a increases, especially for D^. 
The direction of this d-vector is not compatible with the Rashba SO interaction 
at all. For this case, the factor G{k)G{—k) in eq.Q can be negative on a wide 
area of the Fermi surface. Therefore, the direction of the d-vector strongly tends 
to align in the xy-plane. On the other hand, the superconductivity described by 
the gap function D^^ is most stable for large a, because Cq in the Rashba SO 
interaction, which tends to direct the d- vector so that the condition D \\ jCq is satisfied, 
belongs also to Ai irreducible representation. In this sense, D^' is, to some extent, 
compatible with the Rashba SO interaction, though it does not yet fully optimize 
the Rashba interaction, leading to strong inter-band pair correlations. In the case of 
D^^ ~ [(acosA;^. — b cos ky) sin kyX — (bcoskx — a cos ky) sin k,j.y], the parameters a and 
b in the gap function changes as a is increased. For a = 0, {a,b) oc (1,0). When a is 
turned on, a is decreased while b increased so that D^^ would become close to the form 
compatible with Cq. This change is continuous with respect to a and the gap function 
D^^ is transformed gradually. As will be shown later in figure [5] and figure El D^^ has p- 
wave like character for small a and is gradually changed into the /-wave like gap function 
(i.e. a = 6) as a is increased. As seen in figure El Tc for D^' has a minimum around 
a ~ 0.04 and a hump around a ~ 0.06. This a-dependence is understood as follows. As 
mentioned before, is determined by competition and interplay between the pairing 
interaction and the Rashba SO interaction. As a increases, because of the change of 
the electronic structure due to the Rashba SO interaction, the pairing interaction in the 
p-wave channel becomes weak in our model. This results in the overall decrease of Tc for 
D^^ state. On the other hand, the increase of a also changes the structure of the D^^ 
gap function more compatible with the Rashba SO interaction, suppressing inter-band 
pairings between the SO split bands and also associated pair-breaking effects. A slight 
increase for 0.04 < a < 0.06 is caused by this suppression of the inter-band pairings. 
In our model, for a > 0.1, the decrease of Tc for D^^ is substantial. Thus, the pairing 
state with a = b, i.e. /-wave state, can not be realized. 

3.3. Pairing state and transition temperature in the case with parity-mixing 

The pairing states obtained in the previous subsection is drastically changed once we take 
into account the parity mixing of the singlet and the triplet gap functions. According 
to the behaviors of Tcs and in figure [3l the ci-wave and the p-wave pairing states may 
be mixed through the Rashba SO interaction. For the admixture of the gap functions, 
however, only the gap functions which belong to the same irreducible representation of 
the point group are allowed to coexist. The p-wave gap function D^^ with the highest 
belongs to Ai representation of C^^ and the d-wave gap function Dq ~ (cos k^ — cos ky) 
to Bi. This implies that these two gap functions can not be mixed. Then, the next 
candidate is the admixture of the d-wave state and the p-wave state with D^'. The 
symmetry argument allows this admixture. Indeed, we found that the only one solution 
of eq.Q is the d-\-p(Bi) wave state. Figure|l]shows the transition temperature for this 
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Figure 4. a versus Tc for the singlet and triplet mixed Bi symmetric state at U — 5.5. 




Figure 5. The gap functions for the d+p(Qi) state plotted on {k^^ ky) plane. DQ{iujn = 
iirTck) (left panel) and Di{iujn = ii^Tc^k) (right panel) for a — 0.005, C/ = 5.5. 
Green dotted lines indicate locations of the SO split Fermi surfaces. Black dotted lines 
indicate locations of gap nodes where the gap amplitude vanishes. 



d + p{Bi) wave superconductivity as a function of a. For small a, the superconductivity 
is dominated by the spin-triplet state, while for large a it is dominated by the spin- 
singlet state. For intermediate values of a, the two gap functions are strongly mixed 
with the same order of magnitude. 

We show the k dependence of DQlinTc, k) and Di{i7cTc, k) for the d + p(Bi) state 
in figure [5] for a = 0.005 and in figure [6] for a = 0.1. Note that D2{iuJn, kx, ky) = 
Di{iu!n, ky, kx) is satisfied for the Bi symmetric superconducting state. For small a, 
Di exhibits a conventional p-wave behavior, in the sense that the Fermi surfaces cross 
the nodal lines of Di only near (±7r, 0) and the amplitude of Di is much larger than 
that of Dq. Meanwhile, for large a, Di is more like a /-wave state in the sense that the 
Fermi surfaces cross the nodal lines of Di around (±0.47r, iO.Tvr) in addition to (±7r, 0), 
though the single-particle energy is fully gapped with no nodal lines as will be clarified 
in the next section. The change from the conventional p-wave like gap function to the 
/-wave like one is continuous, and actually the /-wave like state should be classified 
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Figure 6. The gap functions for the d+p{Bi) state plotted on {kx, ky) plane. Do{iujn = 
inTc,k) (left panel) and -Di(jw„ = zttTc, fc) (right panel) for a = 0.1, U = 5.5. Green 
dotted lines indicate locations of the SO split Fermi surfaces. Black dotted lines 
indicate locations of gap nodes where the gap amplitude vanishes. 

as a p-wave state with higher harmonics. In our model, even for large a > 0.1, D^^ 
does not change to the form DqjCq (i.e. a genuine /-wave state) for which no inter- 
band pairing exists. This is because Vt does not favor such a structure of the rf-vector. 
The two factors for the determination of the ci-vector, the pairing interaction and the 
Rashba SO interaction, generally have different origins and favor different types of gap 
functions. Therefore, our results imply that, for the case that the pairing interaction 
in a spin-triplet channel is sufficiently strong, it is rather generally hard for the gap 
function to be compatible with the Rashba SO interaction, resulting in the existence of 
the inter-band pairing, in contrast to previous studies on simple models in which it is 
assumed that the spin-triplet pairing interaction compatible with the asymmetric SO 
interaction always exists [311 EHl [36] . In the case that the pairing interaction for the 
triplet superconductivity is very small compared with that for the singlet one and the 
asymmetric SO interaction, the triplet component is induced by the singlet component, 
hence D \\ Cq can be satisfied. 

We note that the relative phase of D to Dq is determined through the Rashba SO 
interaction. The eigenvalue of eq.([9]) for A = (Dq — D • (j)ia2 where (Dq, D) is the gap 
function illustrated in figures [5| and [6] is smaller than that for A = {Dq + D ■ cr)ia2- 
The resulting gap function has no degrees of freedom with respect to the relative phase 
between the singlet component Dq and the triplet component D. 

As mentioned in the previous subsection, the d+p{Qi) pairing state is time-reversal 
invariant. Thus, the asymmetric SO interaction restores the time-reversal symmetry of 
the superconducting state, which is broken in the bulk of Sr2Ru04 where the chiral p+ip 
state is realized. This restored time-reversal symmetry bears an important implication 
for the realization of time-reversal invariant topological superconductivity, as will be 
discussed in the next section. 
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4. A possible realization of topological superconductivity 

In this section, we present the scenario of the Z2 topological superconductivity on the 
basis of the results obtained in the previous section. As mentioned in the introduction, 
the stability of the topological superconductivity is ensured by time-reversal symmetry 
and the existence of a full energy gap which separates the topologically nontrivial ground 
state from excited states [TOl [121 El there should be no nodal lines of the gap, 
from which gapless excitations may emerge, destabilizing the gapless edge modes and 
destroying the topological phase. The time- reversal invariance is evident for our d+p{Qi) 
state, as mentioned in the previous section. Thus, we, here, examine whether the single- 
particle energy for the d + p(Bi) wave state is fully gapped, and there is no nodal lines 
of the gap. To simplify our analysis, we assume the BCS mean field Hamiltonian, which 
is, in the matrix form, 

Ek + aCoik) ■ <7 Afc 

k\ -Sk - aCQ{-k) ■ ct'^ 

with 

Afc = Do{k)i(X2 + D{k) ■ (Tia2, (12) 

where -Do(^) is the spin-singlet gap, and D{k) is the d-vector for the spin-triplet 
component. The energy eigen value of ( ITTl) is obtained as [H] 

Ek± = [el + a'\Co{k)\' + \D{k)\' + Do{kf 

± 2^{ekaCo{k) + Do{k)D{k)Y + {aCo{k) x D{k)Y]^, (13) 

and the minus branch of the eigen values — -Efc±- When D{k) \\ CQ{k), the energy 
spectrum f|T3l) is reduced to Ek± = \Je\^ + A^(fc), with ek± = ^fc ± a\CQ{k)\ and 
A-i-(fc) = Doi^k) ± \D{k)\. The energy spectrum is diagonal with respect to the SO- 
split-band-index ±, and there are no inter-band Cooper pairs. For the p + d wave 
state obtained in the previous section, D{k) || JColk) does not hold for a wide range 
of parameters. In this situation, the energy spectrum (|T3l) is not diagonal with respect 
to the band index, which implies that there exist inter-band Cooper pairs as well as 
intra-band pairs. The condition for the existence of gapless excitations, Ek± = 0, is 
recast in, 

el - a'\Co{k)\' + \D{k)\' - D,{kf = 0, (14) 

ekDo{k)-aC^{k)-D{k) = {]. (15) 

fc-points at which the excitation energy gap closes should satisfy both eqs. (fT4|) and 
(|T5l) . We examine these conditions numerically for the d + p{Qi) state. The calculations 
presented in the previous sections are valid only for T > T^. Thus, for the evaluation of 
f|T^ which requires the magnitude of the gap function, we assume that the maximum 
values of the superconducting gaps Dq and \D{k)\ at T = are obtained from the 
BCS mean field relation A/T^ = 1.764. Using this approximation, we derive fc-points 
satisfying eqs. (fT^ and (fT5i) for the (i + p(Bi) state. The results for some values of a are 
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Figure 7. fc-points satisfying eqs. pi)) (blue line) and (fT5|) (green line) plotted for 
a = 0.005 (left panel) and a = 0.02 (right panel). For a = 0.005, the left-hand side of 
(fT4)) is larger than in the entire Brillouin zone. 



shown in figure [71 We found tliat wfien a is sufficiently small, the left-hand side of 
is positive for all k in the entire Brillouin zone, and thus there are no gapless excitations. 
As a is increased from 0, the left-hand side of f|T4l) decreases, and when a reaches to a 
value ao ~ 0.02, eg. (1141) is fulfilled in a certain region of k where kx ~ i/cy ztO.GSvr. 
The condition (fT5l) is also satisfied exactly on the line = ^ky, because of the Bi 
symmetry; i.e. Di{—ky,kx) = D2{kx,ky) and DQ{kx,kx) = 0. This implies that the 
gap collapses at certain fc-points on the line kx = ky. Our numerical data for the gap 
function indicate that this gap-closing does not occur for < a < Oq. For this parameter 
region, the bulk excitations have the full energy gap in the whole Brillouin zone, which 
is a necessary condition for the realization of the topological superconductivity. 

We, next, consider the adiabatic deformation of the above d -|-p(Bi) wave state to 
a topologically equivalent state. This deformation is achieved by changing parameters 
of the Hamiltonian without closing the bulk gap jTOl HI]. Since any gradual changes 
of parameters cannot change the nonzero topological number which is discrete, the 
existence of the bulk gap ensures the topological stability of the state. As mentioned 
above, for < a < ao, el + \D{k)\'^ > a'^\Co{k)\'^ + Do{ky > holds. Thus we can 
change adiabatically the magnitude of the spin-singlet gap -Do(^) and the strength of 
the SO interaction a to zero without closing the excitation gap. After this deformation, 
the system is equivalent to a combined system oi a p + ip state and a. p — ip state, 
which indeed exhibits the Z2 topological superconductivity [El [ISl [T71 [H] . As a result, 
the d + p{Bi) wave state obtained in the section 3 is topologically equivalent to the Z2 
topological superconductivity. 

In this Z2 topological phase, there are counter-propagating gapless edge states, 
which are Majorana fermions [TH [T8] . The Majorana edge states may give rise to 
intriguing transport phenomena associated with spin currents ^2]. Since there is a 
diamagnetic supercurrent on the boundary surface, for the detection of the spin current 
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carried by edge quasiparticles, thermomagnetic effects may be utilized [IB]. Also, the 
gapless edge quasiparticles may be observed as a zero bias peak of tunneling currents 
[I71[T8]. 

5. Discussion and summary 

In this paper, we have investigated pairing states realized at the (001) interface of 
Sr2Ru04 by using microscopic calculations based on the Kohn-Luttinger-type pairing 
mechanism. It is found that at the (001) interface of Sr2Ru04, the strong admixture 
of p-wave pairings and (i-wave pairings realizes, and thus this system is suitable for 
the exploration of strong parity-mixing of Cooper pairs caused by broken inversion 
symmetry. An important implication of our results is as follows. When there are strong 
spin-triplet pairing correlations in NCSC, the frustration between the pairing interaction 
and the asymmetric SO interaction occurs quite generally, because of incompatibility 
between the pairing interaction and the symmetry of the asymmetric SO interaction. 
This yields substantial spin-triplet inter-band pairs between electrons in two SO split 
bands, even when the size of the SO split is considerably larger than the superconducting 
gap. Because of this feature, the most stable parity-mixed pairing state realized at the 
(001) interface of Sr2Ru04 is the d + p{Qi) wave state, in which time-reversal symmetry 
is restored, in contrast to the bulk Sr2Ru04, which is believed to be in the chiral 
p + ip state with broken time-reversal symmetry. Another intriguing conclusion drawn 
from our results is that this d -|-p(Bi) wave state can be a promising candidate of the 
recently-proposed Z2 topological superconductivity. That is, the d + p(Bi) wave state 
is topologically equivalent to the state that consists of a p + ip state for || pairs and 
a. p — ip state for J,| pairs, which supports the existence of counter-propagating gapless 
edge states carrying spin currents. This feature may be observed experimentally in the 
transport properties for heat currents and spin currents, as discussed in some literature 
[m [TBI 112]. 

Some concluding remarks are in order. For the spin-triplet pairing state obtained 
in the above analysis, inter-band pair correlation between the SO split bands is 
substantially large. This result for inter-band pairings implies that pairing states 
with center-of-mass momentum such as the Fulde-Ferrel-Larkin-Ovchinikov (FFLO) 
state [ISl m] may be more stabilized for some parameter regions compared to uniform 
states considered in the current paper. The stability of the FFLO state depends on 
competition between pairing interaction for this state and the cost of the kinetic energy 
due to the finite center-of-mass momentum. The issue of a possible realization of the 
FFLO state at an interface of Sr2Ru04 is quite intriguing, and should be addressed in 
the near future. 

In the argument for the realization of the topological superconductivity presented 
in this paper, we did not consider effects of the bulk SO interaction due to the c?-electron 
orbitals, but simply assumed that the asymmetric SO interaction overwhelms the bulk 
SO interaction in the vicinity of the interface. Actually, there should exist a domain 



Pairing state at an interface of Sr2RuOi 



16 



boundary between the chiral p+ip state in the bulk governed by the bulk SO interaction 
and the time-reversal-invariant p + d state near the interface. It is highly nontrivial 
how the interaction between these two states affects the stability of the topological 
superconductivity. However, it is expected that as long as the thickness of the region, 
where the asymmetric SO interaction is dominant, is sufficiently large compared to the 
coherence length of Cooper pairs, the topological superconductivity should be stable in 
the vicinity of the interface. Another point required for the realization of the topological 
superconductivity is the fabrication of the interface where the electronic structure is not 
so different from the bulk one. In the (001) surface of Sr2Ru04, however, it is known 
that a structural phase transition occurs and the surface state is ferromagnetic [15] . 
To observe the time-reversal symmetric superconductivity, a very carefully fabricated 
sample without such a structural phase transition is needed. 
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